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Abstract. The rapid neutron capture process (r-process) is thought to be responsible for the creation of more than 
half of all elements beyond iron. The scientific challenges to understanding the origin of the heavy elements beyond 
iron lie in both the uncertainties associated with astrophysical conditions that are needed to allow an r-process to occur 
and a vast lack of knowledge about the properties of nuclei far from stability. There is great global competition to 
access and measure the most exotic nuclei that existing facilities can reach, while simultaneously building new, more 
powerful accelerators to make even more exotic nuclei. This work is an attempt to determine the most crucial nuclear 
masses to measure using an r-process simulation code and several mass models (FRDM, Duflo-Zuker, and HFB-21). 
The most important nuclear masses to measure are determined by the changes in the resulting r-process abundances. 
Nuclei around the closed shells near N=50, 82, and 126 have the largest impact on r-process abundances irrespective of 
the mass models used. 
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Basic properties of nuclei, such as their binding energies 
per nucleon allow the synthesis of the elements up to approx- 
imately iron via fusion reactions in stars from the lightest el- 
ements created by the Big Bang. However, the abundances of 
elements in our solar system contain a substantial number of 
nuclei well beyond iron 1TJ|2][3). The origins of these nuclei are 
entangled in complexity since the heavier elements are thought 
to be made via both slow- and rapid- neutron-capture processes 
(s- and r-processes) [4|. The s-process leads to a network of nu- 
clei near stability while the r-process allows the production of 
nuclei with increasing neutron numbers much further from sta- 
bility, producing neutron-rich nuclei. The astrophysical scenar- 
ios in which the s-process can take place have been identified, 
but a potential site for the r-process is still unresolved|5|. The 
challenge for astrophysical science today is to understand the 
conditions that would provide a major abundance of neutrons 
and lead to successive captures before the nucleus has a chance 
to decay; while on the nuclear side, the challenge is to deter- 
mine the physics of nuclei far from stability where the range 
and impact of the nuclear force is less well known (6][5]]. 

There have been a number of astrophysical scenarios sug- 
gested as possible sites for the r-process. Some of the most 
promising sites include the neutrino driven wind from core- 
collapse supernovae J7J, two-neutron star-mergers |8|, gamma- 
ray bursts [9 1, black-hole neutron star mergers ifTOl . relativistic 
jets associated with failed supernovae [ 1 1 1 or magnetohydro- 
dynamic jets from supernovae fl2l . 

The r-process proceeds via a sequence of neutron captures, 
photodissociations and (3 decays. Simulations of the r-process 
therefore require tabulations of /3-decay lifetimes, neutron cap- 



ture rates and neutron separation energies; photodissociation 
rates are determined from the capture rates and separation en- 
ergies by detailed balance 1131 : 
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In the above expression, T is the temperature, (av) ez,A-x) is 
the thermally-averaged value of the neutron capture cross sec- 
tion for the neighboring nucleus with one less neutron, and 
S n (Z, A) is the neutron separation energy — the difference in 
binding between the nuclei (Z,A) and (Z,A — 1). Nuclear 
masses are crucial inputs in theoretical calculations of each of 
these sets of nuclear data. 

One way to assess the role of nuclear masses in the r-pro- 
cess is to choose two or more mass models, calculate all of 
the relevant nuclear data with the mass model consistently, and 
then run r-process simulations with the different sets of global 
data. Such comparisons are quite valuable and examples in- 
clude Ref. 1 14 15 16 1. Our approach here is quite different. We 
instead focus on the sensitivity of the r-process to the individual 
neutron separation energies within a given mass model, as they 
appear in Eqn.Q] in an attempt to determine the nuclei that have 
the greatest impact on the overall r-process abundances and, in 
turn, identify the most crucial measurements to be made. This 
is the first time that such an attempt has been made and the re- 
sults could potentially be of great significance to both nuclear 
and astrophysical science. 

The study of radioactive nuclei far from stability approach- 
ing the r-process path is one of the global research frontiers 
for nuclear science today. New facilities are being developed 
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Fig. 1. Comparison of the separation energies from Duflo-Zuker |20|, 
HFB-21 1211 . and the experimental masses from 1221 to the FRDM 
1231 values for the tin isotopes. 



in the USA (CARIBU at ANL, NSCL and FRIB at MSU), in 
Europe (ISOLDE at CERN), in France (SPIRAL II at GANIL), 
in Finland (Jyvaskyla), in Germany (FAIR at GSI Darmstadt), 
in Japan (RIKEN), in China (BRIF,CARIF in CIAE Beijing), 
and in Canada (ISAC at TRIUMF). The overarching question 
for this global effort in nuclear science is which measurements 
need to be made IfTTll . 

This study used a fully dynamical r-process nuclear net- 
work code [ 1 8 1 . Inputs to the simulation code include a seed 
nucleus, neutron density, temperature and dynamical timescale 
descriptive of a given astrophysical scenario. In addition, j3 de- 
cay rates, neutron capture rates and neutron separation ener- 
gies are the inputs for the nuclear properties. The simulation 
processes neutron captures, photodissociations, /3-decays, and 
/3-delayed neutron emissions from the start of the r-process 
through freezeout and the subsequent decay toward stability 
fl9l . Fission, while important in some astrophysical scenarios, 
is not significant for the the conditions used here and so is not 
included. 

All the calculations are done for the same initial astrophys- 
ical conditions. The astrophysical scenario used in our simula- 
tions was based on the H or high frequency r-process suggested 
by Qian et al. [24], with an initial temperature of Tg = 1.5 and 
an initial density of 3.4 x 10 2 g/cm 3 . We take the temperature 
and density to decline exponentially as in ll25l with a dynami- 
cal timescale of 0.86 s. While Qian specifies a seed of 90 Se and 
a neutron to seed ratio (N n /N see d) of 86 [24|, here a lighter 
seed of 70 Fe is chosen, which results in N n /N see d = 67 when 
the electron fraction is kept consistent with Qian (Y e = 0.190). 

The nuclear data inputs include beta decay rates from [26 1 
and neutron capture rates from [27|, both calculated with Finite 
Range Droplet Model (FRDM) masses. The measured values 
of S n come from the Audi Mass Evaluation 2003 ll22l . For the 
remaining nuclei, we used the S n values resulting from the cal- 
culated mass values in the FRDM [23]. We subsequently varied 
these theoretical S n for one nucleus at a time by ±25%. In each 
case, the resulting r-process abundance curves were generated 
and compared against the baseline abundances resulting from 
the unchanged S n value. 



Fig. 2. Final r-process abundances for the baseline H-scenario [ 24] 
with 70 Fe seed (black line) compared to simulations in which the neu- 
tron separation energy of 138 Sn is increased (red long-dashed line) 
or decreased (blue short-dashed line) by 25%. The calculated abun- 
dances are normalized to the solar r-process abundances of Sneden et 
al. (4) (points) at A = 130. 

The 25% variation of separation energies was chosen some- 
what arbitrarily. A comparison of the ratio of separation en- 
ergies extracted from measured masses or theoretically calcu- 
lated separation energies with the FRDM calculated values is 
shown for the Sn isotopes in Figure Q] This indicates that the 
25% value is a reasonable variation estimate far from stability. 

An example of the resulting abundance patterns is shown 
in Figure |2 where the baseline pattern is compared to the final 
abundance patterns produced by simulations in which the sep- 
aration energy of 138 Sn was increased or decreased by 25%. 
This comparison can be quantified by summing the differences 
in the final mass fractions: 

F± = l00j2\ X ba S eline{A)-X ±AS „(A))\, (2) 
A 

where X (A) = AY (A) is the mass fraction of nuclei with 
mass number A (such that J^a X(A) = 1), and the sum of A 
ranges over the entire abundance curve. This quantity is largest 
when the curves differ near the peak abundances, giving pref- 
erence to those regions. 

The values of F = (F+ + F-)/2 are calculated for 3010 
nuclei from 58 Fe to 294 Fm. Figure [3] shows the nuclei whose 
separation energy variations result in the greatest changes in the 
resulting r-process abundances. Nuclei that have the greatest 
impact on the r-process are those neutron rich nuclei near the 
closed shells at Z = 28 and 50, and N = 50, 82, and 126. 

A natural question to ask is the dependence of these results 
on the mass model used. Therefore, similar calculations were 
performed using four additional mass models, the Duflo-Zuker 
(DZ) [20], the Extended Thomas Fermi plus Strutinsky Inte- 
gral with shell Quenching (ETFSIQ) [30], the Hartree-Fock- 
Bogoliubov (HFB-21) |[2D . and the F-spin [31] model in addi- 
tion to the FRDM. All models take advantage of very different 
physics ingredients to calculate the masses of nuclei far from 
stability. Each of the calculations performed started with the 
same initial astrophysical conditions and again varying individ- 
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Fig. 3. Comparison of the sensitivity to mass values determined by 
Equation 2. The separation energies far from stability were generated 
by the FRDM (23), Duflo-Zuker (20), and HFB-21 (21). The scale is 
from white to dark red, indicating regions with a small change to a 
substantial change in the resulting abundances. For reference, stable 
nuclei have been included as black crosses and the magic numbers 
have been indicated by thin lines. Superimposed on the sensitivity re- 
sults are the limits of accessibility by CARIBU 1 28 1 and the proposed 
FRIB intensities |29|. In both cases, we have plotted the conservative 
limits of what can be produced and measured in mass measurements. 



ual separation energies by ±25%. The results are astounding. 
In each case, the nuclei with the greatest impact were generally 
the ones near the major closed shells independent of the cho- 
sen mass models. Figure[3]shows the resulting sensitivity plots 
from three of the mass models; the FRDM, DZ, and HFB-21 
models. Nuclei near the closed shells of N=50, 82, and 126 rise 
above all the others in impact. The nuclei with the most impact 
on the r-process abundances cluster around 132 Cd and 138 Sn. 
In this region, the nuclei are 131 " 134 Cd, 132 " 137 In, 135 " 140 Sn, 
139 ' 141 Sb. There are also specific low mass nuclei such as 82 Cu, 
85 Zn, and 88 Zn that are important. 

In trying to understand these results, we know that there 
are two ways that an individual neutron separation energy can 
influence the r-process abundance distribution. The first is a 
long-recognized 1132) equilibrium effect, and the second is an 
early-freezeout photodissociation effect, recently pointed out 
in 1133) . In the classic view, the r-process takes place in condi- 
tions of (n, 7)-(7, n) equilibrium, where abundances along an 
isotopic chain are determined by a Saha equation: 
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Fig. 4. Shows the mass fractions of 136 Sn (purple), 138 Sn (blue), and 
140 Sn (aqua) for the baseline r-process simulation (top panel) and 
the simulation with the separation energy of 138 Sn decreased by 25% 
(middle panel). The bottom panel compares the neutron abundance for 
the two simulations (black and red lines, respectively). 
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where the Gs are the partition functions, N n is the neutron 
number density, and m„ is the nucleon mass. The relative abun- 
dances of the different isotopic chains are then determined by 
the /3-decay lifetimes of the most populated nuclei along each 
chain. As described in Eqn.[3] any change to an individual sep- 
aration energy will cause a shift in the abundances along the 
isotopic chain. This can have a global impact on the final abun- 
dance pattern, particularly if the affected nucleus is highly pop- 
ulated and material is shifted to a nucleus with a significantly 
faster or slower /3-decay lifetime. For example, consider the 
case of 138 Sn, a nucleus just above the N = 82 closed shell re- 
gion. In the baseline simulation, 138 Sn is the most abundant tin 
isotope, and 136 Sn 140 Sn are much less abundant. Their mass 
fractions are shown as a function of time in Fig.[4ja); their rela- 
tive values follow those predicted by Eqn.[3]until about t ~ 1 .2 
s, when equilibrium begins to fail and the nuclei primarily /3- 
decay to stability. If the simulation is repeated with neutron 
separation energy of 138 Sn reduced by 25%, we see that the 
equilibrium abundance of this nucleus is drastically reduced, 
as expected from Eqn. [3] and shown in Fig. Sfb). Material is 
instead shifted to 136 Sn, which has a /3-decay lifetime approx- 
imately 1.6 times that of 138 Sn (and 5.3 times the lifetime of 
140 Sn, which is also depleted by the shift). As a result, more 
material is stuck in the tin isotopic chain compared to the base- 
line simulation, and the overall rate at which neutrons are con- 
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sumed is slowed, as shown in Fig.|4jc). This impacts the avail- 
ability of neutrons for the whole abundance pattern and results 
in changes throughout the pattern. The second mechanism, in 
contrast, operates once (n, r y)-( r y, n) equilibrium begins to fail, 
and individual neutron capture and photodissociation rates be- 
come important. Since the neutron separation energy appears in 
the exponential in Eqn.Q] photodissociation rates are quite sen- 
sitive to this quantity. Changes in individual photodissociation 
rates during freezeout can produce local shifts in abundances, 
which can translate into global abundance changes if they alter 
the late-time availability of free neutrons. This mechanism is 
described carefully in ll33l . Odd-TV nuclei, which tend to be in 
equilibrium only briefly if at all, are particularly susceptible to 
these non-equilibrium effects. 

In conclusion, this study of 3010 nuclei via an r-process 
simulation tested the sensitivity of the r-process abundance yie- 
lds to the theoretical mass values of neutron rich nuclei present- 
ly unknown in the laboratory from several different mass mod- 
els, the results are shown here for three of them (FRDM|23|, 
Duflo-Zuker ll20l . and HFB-2 H12D ). The results are uniform 
and conclusive in highlighting the importance of nuclei near 
closed shells. Essentially the same set of nuclei emerge as hav- 
ing the highest impact on the r-process irrespective of the vary- 
ing physics ingredients of the different mass models. The nu- 
clei with greatest impact on the r-process — neutron rich iso- 
topes of cadmium, indium, tin, and antimony in the N = 82 
region, nickel, copper, zinc, and gallium in the N = 50 region, 
and thulium, ytterbium, lutetium, and hafnium in the N = 126 
region — should be of highest priority to measure in the vari- 
ous exotic beam facilities around the world. Table Q] shows the 
top 25 nuclei with the greatest impact on the r-process for the 
three models. Since the particular isotopes of these elements 
that have the greatest impact can shift depending on the astro- 
physical conditions, a future paper will explore the effects of 
various astrophysical scenarios in determining the most impor- 
tant nuclei to measure. 

This work was supported by the National Science Founda- 
tion through grant number PHY0758 100 and the Joint Institute 
for Nuclear Astrophysics grant number PHY0822648. 
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